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ABSTRACT 

The inspiral of a turbulent molecular cloud in the Galactic centre may result in the 
formation of a small, dense and moderately eccentric gas disc around the supermassive 
black hole (SMBH). Such a disc is unstable to fragmentation and may lead to the 
formation of young massive stars in the central parsec of the Galaxy. Here we perform 
high-accuracy direct summation TV-body simulations of a ring of massive stars (with 
initial semi-major axes 0.1 ^ a/pc ^ 0.4 and eccentricities 0.2 ^ e ^ 0.4), subject 
to the potential of the SMBH, a stellar cusp, and the parent gas disc, to study how 
the orbital elements of the ring evolve in time. The initial conditions for the stellar 
ring are drawn from the results of previous simulations of molecular cloud infall and 
disruption in the SMBH potential. While semi-major axes do not evolve significantly, 
the distribution of eccentricities spreads out very fast (s=s 1 Myr) as a consequence 
of cusp precession. In particular, stellar orbits with initial eccentricity e > 0.3 (e < 
0.3) tend to become even more (less) eccentric, resulting in a bimodal eccentricity 
distribution. The distribution is qualitatively consistent with that of the massive stars 
observed in the Galactic centre's clockwise disc. 

Key words: black hole physics - methods: N-body simulations - stars: kinematics 
and dynamics - Galaxy: centre 
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1 INTRODUCTION 

Despite the strong tidal field exerted by the super massive 
black hole (SMBH) in the Galactic centre (GC, iMorrisI 
1993), more than a hundred young massive stars are ob- 
served in the central parsec of the Galaxy. 

The majority are O-type and Wolf-Rayet (WR) stars, 
of which about two thirds show clockwise motion when 
projected on the plane of the sk y while the remain- 
der show countercloc k wise motion (IPaumard et al. 2006; 
Navakshin et alJIgOOfj ; iMartins et all 120071 : iLu et allhoOS 



Bartko et all 120091 '). The WR/O stars in the larger group 



orbit the SMBH in a thin disc, usually referred to as the 
clockwise (CW) disc, with an estimated age of ~ 6 Myr 
and an average eccentricity of ~ 0.4. The smaller group 
of stars on counter-rotating orbits ma y indicate the pres- 
ence of an additional dissolv ing disc l|Bartko et al.l 120091 ; 
lLockmann fc Baumg ardt 2009]) . The CW disc shows a sharp 
inner truncation at about one arcsecond (or 0.04 pc). 
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The S-cluster iSchodel et al.l 120031: iGhez et al.l 12001 



120051 : lEisenhauer et al.ll2005l ; iGillessen et al.l 120091 1 located 
in the innermost arcsecond is a group of about 20 B-type 
stars moving on eccentric and randomly oriented orbits, with 
estimated ages of 20 — 100 Myr. The origin of the S-cluster 
is largely controversial, with main models including the in- 
spira l of a young star cluster ([Gerhard 200lj^ Kim fc Morris 



20031 : iMcMillan fc Portegies Zwartl 12001 iGiirkan fc Rasic 
20051 : iLevin et al.ll2005l : iFuiii et all 1120081). possibly hosting 
an in te rmediate-mass black hole ((Hansen & Milos avlieviq 
20031: iBerukoff fc Hansen] 120061: iMerritt et al.l [20091 ; 

20101 : 



Gualandris fc MerrittJ 
'20091 1 



20091 : iGualandris et al 



Fuiii et al 
on l ow- angular 



and the tidal disrupti on of binaries 
momentum orbits (iGould fc Quillenl 



20031: IPerets et al.l [20071: iPeretsI [20091: IPerets et all 120091 : 



Antonini et al.ll2010l : iPerets fc Gualandrid BoiOl 



The formation of the CW disc, which represents the 
target of this work, is also debated, but is often at- 
tributed to in-situ star formation from an infalling gas cloud 
jMorris! Il993l: iNavakshin fc Cuadral 120051: [Bonnell fc Ricel 



120081 : iMapelli et al.l 120081 : IWardle fc Yusef-Zadehl 120081 : 
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Hobbs fc Navakshinl 120091; Alig ct al. 1 l201ll ; iMapelli et alj 
20121 ). In particular. iMapelli et al.1 (|2012h . hereafter M12, 



perform A?-body/Smooth Particle Hydrodynamics (SPH) 
simulations of the infall of a turbulent molecular cloud to- 
ward the GC and study the formation of stars in the region. 
They find that the tidal disruption of the cloud results in 
the formation of a small (radius < 0.5 pc), dense (> 10 s 
atoms cm~ 3 ) and eccentric (e ~ — 0.6) gas disc around 
the SMBH. The simulations by M12 are the first following 
the fragmentation of the gas disc into self-bound clumps. 
Stars originate from such clumps in a ring at a distance of 
0.1 — 0.4 pc with moderately eccentric orbits (e ~ 0.2 — 0.4). 
The mass function of the stars is top-heavy if the local back- 
ground temperature is sufficiently high 100 K) and if 
the parent cloud is sufficiently massive (> 1O 5 M ). These 
propertie s are in good agreement with observati ons of the 
CW disc (|Paumard et al.l2006l ; lBartko et al-hoOSh . and lend 
further support to the in-situ formation model for the disc. 
The simulations by M12, however, follow the evolution of 
the system for only about 0.5 Myr since the formation of 
the gaseous disc, whereas the estimated age of the CW disc 
is of about 6 Myr. 

Here we use high-accuracy direct summation iV-body 
simulations to follow the evolution of the stellar ring for 
10 Myr under the effect of the combined gravitational po- 
tential due to the SMBH, the stellar cusp, and the residual 
gas disc. We stress that this is the first calculation ever per- 
formed to follow both the short term star formation phase 
and the long term evolution of the disc in a consistent man- 
ner. We find that the disc retains its initial inner and outer 
radius, and that stars remain largely confined to the ini- 
tial ring plane, which is aligned with the gas disc. Inter- 
estingly, we obs erve evidence of the e ccentric disc instabil- 
ity described by iMadigan et"aH l|2009h . according to which 
precession of the orbits due to the presence of the stellar 
cusp induces coherent torques that drive orbital eccentric- 
ities away from their initial values. As a result, a bimodal 
eccentricity distribution is established which is consistent 
with observations of the CW disc in the GC (|Bartko et al.l 
I2009D . 



2 N-BODY SIMULATIONS 

We use the direc t summation iV-body code <^GRAPE 
l|Harfst et al.ll2007f ) modified to include the effects of an ex- 
ternal potential due to the stellar cusp and the gas disc. 
We adopt a mass of M. = 3.5 x 10 6 M B for the SM BH and 
the broken power- law profile of iGenzel etafl (|2003h for the 
stellar cusp density 



p(r) = po(r /r)i 



(1) 



with r = 10" ~ 0.39 pc, p = 1.2 x lO 6 M pc~ 3 , 7 = 1.4 
for r < ro and 7 = 2.0 for r ^ ro. For the residual gas disc 
we adopt a projected density profile of the type 



E(R) = E ( J R /fi)° 



(2) 



with parameters a = 1.14, Ro = 1 pc and £0 = 2.64 x 
10 3 Mq pc -2 derived from fitting a power-law profile to the 
gas remaining after star formation in the simulation of M12. 
The mass profile of the stellar cusp and gas disc are shown in 
Fig-Ill The gravitational potential generated by a thin disc 
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Figure 1. Solid line: mass profile of the stellar cusp whose gravi- 
tational effect is included in the simulations in the form of an an- 
alytic external potential. Points: mass profile of the gaseous disc 
obtained by the SPH simulations of M12. Dashed line: power- 
law fit to the gas disc profile adopted to generate the external 
potential by means of Eq.|3] Dotted line: mass profile of stellar 
ring as obtained by M12 and used as initial conditions for the 
integrations. 



of kn own E(-R) is given by the integral (|Binnev fc Tremaind 
1 19871 , Eq.2-142b) 



<S>(R,z) 
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where 



iRR' 



(R + R'Y + z 2 



(3) 



(4) 



and K(k) represents the complete elliptic integral of the first 
kind. Differentiating Eq.[3]with respect to R and z gives 
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(6) 



where £(k) represents the complete elliptic integral of the 
second kind. The equations for the gravitational potential 
and forces can be solved numerically on a cylindrical coor- 
dinate grid. 

Relevant timescales for the stars are shown in Fig. [2] 
as a function of distance from the SMBH. At distances of 
0.1 — 0.4pc, precession due to general relativity (GR) is not 
important on timescales of 10 Myr, and we therefore neglect 
it in the simulations. On the other hand, Newtonian preces- 
sion due to the stellar mass distributed within the orbits of 
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Figure 2. Relevant timescales for the ring stars in the GC: orbital 
period (solid line), timescale of GR precession (dashed lines) and 
timescale for Newtonian precession (dot-dashed line). The shaded 
area indicates the approximate location of the CW disc. An or- 
bital eccentricity of 0.3 is assumed. 




the stars has a timescale of (e.g. iMerritt et al.|[201ll ) 

M. 



Tnp = 



Ale- 



-T OT b/(e) : 



(7) 



where T or b is the orbital period of a star with semi-major 
axis a, M cusp is the total stellar mass enclosed within the 
orbit and 



/(e) = 



Tnp is of about 10 s yr for the stars in the ring. For 7 = 3/2, 
Tnp is independent of a. For our adopted broken power- 
law profile, Tnp is only weakly dependent on distance (see 
Fig.®. 



2.1 Initial conditions 

We adopt the outcomes of run E by M12 as initial condi- 
tions for the simulations presented in this paper. In particu- 
lar, M12 performed N-body/SPH simulations of a gas cloud 
evolving in the potential of a M. — 3.5 x 10 6 M SMBH and 
a stellar cusp as described in Eq.Q] The impact parameter of 
the cloud with respect to the SMBH is 10 -2 pc and the ini- 
tial velocity is close to the escape velocity from the SMBH 
at the initial distance (25 pc). The cloud is quite massive 
(1.3 x 1O 5 M0), initially spherical (with a radius of 15pc), 
marginally self-bound and seeded with supersonic turbulent 
velocities. The thermodynamic treatment in run E in cludes 
radia tive cooling with accurate recipes for opacity (|Bolevl 
120091 ) and an irradiation temperature Ti rr = 100 K, mim- 
icking the background temperature of the GC. Due to the 
SMBH shear, the molecular cloud is disrupted in « 10 5 yr 
and, because of its low orbital angular momentum, part of 
the gas settles in a very small and dense disc surrounding 
the SMBH (see Fig. [3]). A fraction of the gas in the disc frag- 
ments forming self-bound clumps (with density ) 2 x 10 12 
cm" 3 ). As a combined effect of the initial cloud mass and 
of the floor temperature, the mass function of the clumps is 
top-heavy, ranging between 1 and 60 Mq after 4.8 x 10 s yr. 



Figure 3. Density map of the gas in run E of M12 at t = 
4.8 X 10 5 yr. The simulation was projected in the plane where 
the gaseous disc (at the centre) is seen face-on. The box measures 
20 pc per edge. 



Here we assume that each self-bound clump formed in 
run E of M12 becomes a star, without accreting any further 
gas particles after t — 4.8 x 10 J yfl Thus, we replace each gas 
clump with a single star particle having mass equal to the 
total mass of the clump itself, and position and velocity cor- 
responding to those of the centre-of-mass of the clump. The 
resulting stellar ring, with a total mass of ~ 4.3 x 10 3 Mq, is 
shown in the top panels of Fig. [4] The edge-on view shows 
a small degree of warping. 

In this paper, we cannot integrate the star particles to- 
gether with the gas particles, as it would be prohibitive to 
follow the gas for £ 1 Myr. Thus, we substitute the remain- 
ing gas particles with an analytical potential (Eq. [2} very 
close to the final distribution of gas in run E of M12 (see 
their figure 10). The total mass of the gas disc within 0.4 pc 
(i.e. the outer edge of the stellar ring) is ~ 1O 4 M0. For 
simplicity, the potential is assumed constant in time. While 
this assumption may seem unrealistic, it is justified by the 
fact that the evolution of the ring stars is dominated by the 
SMBH and the stellar cusp, and the gas disc potential repre- 
sents a small correction to the total potential. We performed 
an additional simulation excluding the gas disc potential and 
we found a very similar evolution for the ring stars. 

The SMBH is treated as a point-mass A"-body particle, 
similarly to the ring stars. However, given the small number 
of particles in the simulations, the SMBH would be subject 
to an unr ealistically lar ge motion at the center of the poten- 
tial well <|MerrittJl200"5T ). which might affect the dynamics of 
the stars. For this reason, at each integration step we force 
the SMBH to remain fixed at the center of the reference 
frame. 



1 This assumption is quite conservative for the masses of stars, 
as it is realistic to expect further gas accretion in run E of M12. 
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Figure 4. Snapshots of the ring stars at different times during 
the evolution (from top to bottom t = 0,3, 6 Myr, in a reference 
frame in which the stars lie in the x — y plane. The star symbol 
marks the position of the SMBH. 



3 RESULTS 

We let the ring of stars evolve for 10 Myr under the effect of 
the gravitational potential generated by the SMBH, the stel- 
lar cusp, and the gas disc. Due to the presence of extended 
mass (stars and gas) within the orbit of each ring star, stars 
precess in their orbital plane. Because of its much higher 
mass, the stellar cusp is responsible for most of the preces- 
sion. While the dependence of the precession time on dis- 
tance is rather weak (see Fig. [5} , the dependence on eccen- 
tricity may be significant. This implies a different precession 
time for stars with different eccentricity, and it results in a 
quick loss of coherence in the ring, as can be seen from the 
snapshots at 3 Myr and 6 Myr shown in Fig. [4] Only in the 
absence of a stellar cusp would the ring evolve coherently. 

The orbital elements of all stars are shown in Fig. at 
different times during the integration. While the semi-major 
axes distribution remains unchanged during the evolution, 
the distribution of eccentricities evolves toward a bimodal 
distribution, as can be seen in Fig. [6] While the initial dis- 
tribution of eccentricities shows a single peak at e ~ 0.3, 
the distribution after 6 Myr has a primary peak at e ~ 0.1 
and a secondary peak at e ~ 0.5, with a tail that extends to 
about e ~ 0.7. The cumulative distributions of eccentricities 
are shown in the insert. 

We attribute this behaviour to the eccentric disc insta- 
bility which was pre dicted to occur in ecce ntric discs subject 
to cusp precession jMadigan et all 120091 ). Due to the fact 
that /(e) in Eq. [7] is a monotonically increasing function 
of eccentricity, an orbit with initial eccentricity in excess of 
the average will have a longer precession time, which will 
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Figure 5. Orbital elements of the ring stars at the start of the in- 
tegration (left panel), after 3 Myr (middle panel) and after 6 Myr 
of evolution (right panel). 




Figure 6. Eccentricity distribution of the ring stars at the start of 
the integration (filled area) and after 6 Myr of evolution (hatched 
area). The insert in the top-right corner shows the cumulative 
distribution of eccentricities at different times: at the beginning 
of the integration (dashed line), after 3 Myr (dotted line) and after 
6 Myr (solid line). 



result in a stronger torque from the other stars in the ring. 
The torque will reduce its angular momentum and therefore 
cause an increase in eccentricity. The opposite effect will 
happen to stars with orbits initially more circular than the 
average. The instability is illustrated in Fig. [7] which shows 
the time evolution of the orbital eccentricity for a random 
subset of stars in the ring. While stars initially more eccen- 
tric than the average (~ 0.3) tend to become more eccen- 
tric, stars initially less eccentric than the average tend to 
circularise. The observed timescale for the process is of only 
r^ 1 Myr. This is in a greement with the predicted timescale 
jMadigan et alj r2009) 



Tnp 



M ■ 



1/2 r 



de 



1 

W) 



-,-1/2 



(8) 

where M cusp and M r i ng are the cusp mass and the ring mass, 
respectively. In our case, 30 M cusp / A/ r ing ;$ 60 for 0.1 ^ 
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Figure 7. Eccentric orbit instability in a random subset of stars. 
The timescale for the process is about 1 Myr. 



r/pc ^ 0.4 and T occ ~ 5 — 10 for a circular orbit. As the 
ring loses coherence after this time, torques become weaker 
and changes in the eccentricities become less pronounced. 

This distribution is consistent wit h, though not identi- 
cal to, the one found for the CW disc. iBartko et all (j.2009) 
perform a statistical analysis of the candidate disc members, 
assuming Keplerian orbits in the gravitational potential of 
the SMBH and the stellar cusp, and derive distributions 
for the orbital elements of the stars. While these distribu- 
tions are subject to many uncertainties, the main properties 
of the eccentricity distribution appear robust. In particu- 
lar, they report a bimodal distribution for the reconstructed 
eccentricities with peaks at e ~ 0.3 and 0.9, though they 
caution that the high eccentricity stars may not dynami- 
cally belong to the disc. Once corrected for projection ef- 
fects, the mean eccentricity in the CW disc is ~ 0.37, only 
slightly higher than what we find. The difference is likely due 
to the eccentricity with which the stars form from the gas 
disc, which in turn depends on the orbital parameters of the 
cloud. An average eccentricity of it 0.4 could be produced 
starting with an even more radial orbit for the molecular 
cloud, which might result from cloud-c loud collisions near 
the GC (e.g. iHobbs fe Navakshinl 120091) . Very large eccen- 
tricities ( £ 0.99) like those found bv lMadigan et al.l |2009h . 
however, seem hard to produce even starting from a more ec- 
centric gas disc. This has important implications for models 
of the origin of the S-cluster in the GC. For tidal disrup- 
tion and capture to occur, a binary must come at least a s 
close to the SMBH as the tidal radius. iMadigan etaH l|2009h 
find that for a 4 x 10 3 Mq disc with initial eccentricity of 
0.6 the maximum eccentricity attained by the stars is about 
0.9, and yet there are no stars with a pericentre distance 
small enough to be tidally disrupted. Only discs more mas- 
sive than 10 4 M© and with initial eccentricities in excess of 
0.6 result in a significant number of tidal disruptions. We 
conclude that the S-stars cannot have formed from a disc 
like the one produced here. 

In addition to finding smaller eccentricities than in some 
of Madigan's simulations, we find smaller inclinations with 
respect to the initial ring plane. Inclinations after 6 Myr are 
generally smaller than ~ 10 degrees and the maximum we 
record is ~ 16 degrees. Inclinations do not seem to correlate 



Figure 8. Inclinations of the ring stars after 6 Myr of evolution as 
a function of eccentricity (left panel) and semi-major axis (right 
panel). The inclination is calculated as the angle between the 
angular momentum of the star and the z-axis, which is parallel 
to the initial total angular momentum of the ring. 



with eccentricities (see Fig. [8jl, though we note that if we 
divide the sample of stars in two groups, having e < 0.3 and 
e ^ 0.3, the inclinations in the two groups evolve in a differ- 
ent way. A 2D Kolmogorov-Smirnov test gives probability 
values of 0.07 at the start of the simulation, 0.01 after 3 Myr 
and 2.1 x 10 -6 after 6 Myr that the inclinations in the two 
groups originate from the same distribution. This difference 
is likely due to differential heating between more eccentric 
and less eccentric stars. On the other hand, Fig. [8] shows a 
trend for higher inclinations to develop for stars closer to 
the SMBH. This is a result of the initial warping present in 
the stellar ring, which in turn reflects the warping of the gas 
disc. Larger inclinations might originate from more strongly 
warped gas discs, but additional SPH simulations would be 
required to test this scenario. The rms inclination grows 
slowly with time, roughly as t 1 ^, while the rms eccentricity 
remains approximately constant at e ~ 0.3. This scaling is 
som ewhat natter than w hat is found in lCuadra et all (2008) 
and lHaas fc Subi] (|2012l ), and it may be attributed to the 
different initial conditions adopted in this work, in particu- 
lar to the non zero mean initial eccentricity and/or warping 
of the stellar ring. The modest inclinations achieved by the 
simulated stars imply that the single cloud infall scenario 
cannot account for the high inclination of some of the young 
stars observed in the GC (those claimed to belong to a coun- 
terclockwise disc). A cloud-cloud collision has been shown to 
produce significantly warped stellar discs and additional fil- 
ame ntary structures at large inc linations with respect to the 
disc (|Hobbs fc Navak shin 200|]). Alternatively, the interac- 
tion between two approximately coeval and highly inclined 
stellar discs may provide an explanation for the existence of 
stars at high inclination (|Lockmann fc Baumgard"t]|2009l ). 

Finally, we find a very modest heating of the stellar 
ring during the evolution. This is likely due to the absence 
of heating from the stellar cusp in our simulations, which 
would result in an increase in the mean ring height (see e.g. 
iPerets et aIll200Sj ). 
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4 DISCUSSION AND CONCLUSIONS 

We studied the evolution of the ring of stars formed in the 
GC from fragmentation of the gas disc deposited by an in- 
spiralling molecular cloud. Our calculation is the first of its 
kind to follow, in a consistent manner, both the star forma- 
tion phase in the gaseous disc and the long term evolution 
of the stellar disc. We find that, while the ring retains the 
original distribution of semi-major axes, and therefore also 
the initial inner and outer radius, the distribution of eccen- 
tricities evolves in time due to the onset of the eccentric disc 
instability. Torques exerted by other stars in the ring result 
in a change in the magnitude of the angular momentum and, 
as a consequence, in the eccentricity. As stars evolve away 
from the average eccentricity, a bimodal distribution is es- 
tablished, with a primary peak at e ~ 0.1, a secondary peak 
at e ~ 0.5, and a tail that extends to e ~ 0.7. This is qual- 
itatively consistent with the distribution found for the CW 
disc stars. We predict that a quantitative agreement would 
require a molecular cloud with an initial orbit more radial 
than in run E of M12, but a new suite of hydrodynamical 
simulations is necessary to test this prediction. On the other 
hand, even a better description of the cloud thermodynam- 
ics and of the gas shocks taking place during the infall of 
the cloud might affect the formation and thus the orbital 
pr operties of the young ste llar ring (see e.g. the discussion 
in I Jiang fe Goodmanll201ll . and references therein). 

Our simulations support the scenario in which the thin 
CW disc originated from fragmentation of a moderately ec- 
centric gas disc resulting from the tidal disruption of a tur- 
bulent molecular cloud on a low angular momentum orbit. 
The origin of the WR/O stars at large inclinations, as well as 
those belonging to the S-cluster, however, seems to require 
additional physical processes, which are currently under in- 
vestigation. In the case of the S-stars, relaxation processes 
against the background spherical cusp of stars and remnants 
are unable to explain the observed dis tribution of eccen - 
tricities within the lifetime of the stars |Perets et al.l I2009T ) 
and different mechanisms are required (for a review s ee e.g . 
iPerets fc Gualandris! [2O10l l . Similarly. ICuadra et al.l (|2008T ) 
have shown that relaxation in the disc does not produce ec- 
centricities and inclinations as large as those observed for the 
WR/O stars outside the CW disc. While some outliers may 
be the result of vecto r resonant relaxation betw een the disc 
and the stellar cusp (|Kocsis fc Tremainell20Tll ). several au- 
thors have investigated the possibility of an additional per- 
turbing pot ential, which would pro duce Kozai-Lidov type 
oscillations fKozai 1962; Lidov 196^) in the eccentricity and 
inclination of the stars. In particular. ISubr et al.l (|2009T ) sug- 
gest that precession due to the gravity of the circumnuclear 
disc (CND) at the edge of the SMBH's sphere of influence 
might have driven some of the stars away from the CW 
disc. An axisymmetric perturbation due the presence of an 
additional gas/stellar disc, an intermediate-mass black hole 
or the CND mig ht have observabl e effec ts also on the stars 
in the CW disc. I Lockmann et al.l |2008T l studie d the effect 
of a s econd highly inc l ined stellar disc, while Sub r et ah! 
(2009) and lHaas et al.l 1120111 ) considered the effect of the 
CND. [Ldckmann et al. (2008) found that, in the absence of 
a stellar cluster, precession and Kozai oscillations due to the 
presence of a second, inclined disc result in stars achieving 
high eccentricities. The inclusion of a stellar cusp, however, 



has been shown to damp Kozai osc illations in the disc (e.g. 
I Chanel [20091 ; [Ldckmann et al.ll2009T l. at least in the case of 
an analytic spherical cusp. In this case, the overall influence 
of the cusp can be characterised by a decrease in the ampli- 
tude of the eccentricity and inclination oscillations, and by 
a sho r tening of their period (|Karas fc Subrll2007l ; ISubr et al.l 
20091. IChanel [|2009D examined the effectiveness of the Kozai 
mechanism induced by stellar discs in the Galactic centre 
and found that Kozai oscillations are suppressed if the pe- 
riod for apsidal precession induced by the cusp is shorter 
than the period for Kozai oscillations. In the case of the two 
disc scenario, the mass of the stellar cusp in the GC is over 
an order of magnitude larger than the critical mass required 
to suppress the Kozai mechanism. 

In general, the damping of the Kozai mechanism de- 
pends on the parameters of the perturbing potential and on 
the orbital parameters of the stars, and oscillations might 
be important for certain sets of par ameters. The per iod of 
Kozai oscillations can be written as (|Subr et al.l l2009) 



Tk = 



M. 



Ftp 



M P aVGM.a 



M. ( R P 

Mp 



(9) 



where Mp and Rp are the mass and radius of the perturbing 
axisymmetric potential. If the perturber is a second disc at 
approximately the same distance as the CW disc from the 
SMBH, R P ~ a and T K /T NP w M cusp /M P . If the perturber 
is the CND, then Rp ~ (5 — 20)a. In both cases, it seems 
unlikely that Kozai oscillations might dominate the evolu- 
tion of the CW disc in the presence of a spherical cusp. The 
timescale for Kozai oscillations is also typically longer than 
the timescale for the eccentric disc instability (see Eq.|8]), 
suggesting that Kozai oscillations should not significantly 
affect the development of the eccentricity growth studied 
here, even if this is theoretically possible for certain combi- 
nati ons of parameter s . 

ILockmann et al.1 |2009t) explored the effect of a granular 
JV-body cusp as opposed to an analytic smooth cusp. They 
found that both a smooth and a granular cusp cons istent 
with the one derived for the GC (|Schodel et al.1 120071 ) sup- 
press Kozai cycles in the two disc scenario. However, in the 
case of an A r -body cusp, two-body and resonant relaxation 
effects may temporarily drive the eccentricty of some stars 
to large values. This effect is a function of the cusp gran- 
ularity, i.e. of the mean and maximum stellar masses, and 
is most efficient in the case of a cusp of stellar rem nants. 
Similar results were obtained by lHaas fc Subrl (|2012T) . 

It is not yet clear whether secular evolution can play 
an important role under realistic conditions, i.e. a realistic 
gravitational potential and granularity of the stellar cusp. 
Shape and parameters of the potential of an additional disc 
and/or the CND are not well determined, and the exact 
composition of the GC cusp is not well known. Further in- 
vestigation of these issues is deemed necessary and will be 
presented in an upcoming work. 
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